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PREFACE 


This report details the test direction planning problems associated with 
cryogenic wind tunnels, analyzed as a part of the project "Modeling and Control 
of Transonic Cryogenic Tunnels," sponsored by NASA/Langley Research Center 
(LaRC) under research grant NSG-1503. The report is concerned with realizing 
desired flow Reynolds number-Mach number combinations at which data is sought, 
with minimum liquid nitrogen consumption. The contents of this document 
complement the reports on modeling phase activity in reference 1 and the 
control analysis phase activity in reference 2. 
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MINIMUM ENERGY TEST DIRECTION DESIGN IN 
THE CONTROL OF CRYOGENIC WIND TUNNELS 


By 

S. Balakrishna 1 
INTRODUCTION 

Need for realization of full-flight Reynolds number flows past test 
models in relatively small wind tunnels was keenly felt in the 1970' s. The 
advent of the cryogenic wind tunnel concept is directly attributable to 
this need for high Reynolds number flow in wind tunnels (ref. 3) . The 
cryogenic wind tunnel concept consists of operating the test medium of a 
conventional tunnel at cryogenic temperatures down to 80 K. Nitrogen gas, 
cooled by injected liquid nitrogen, has proven to be an ideal candidate 
for the cryogenic tunnel test medium because of its near perfect behavior in 
insentropic flow (ref. 4). Further, cryogenic operation of a wind tunnel 
results in reduced fan power consumption and no penalty in flow dynamic 
pressure. 

In a cryogenic tunnel, the flow parameters — Reynolds number, Mach number 
and flow dynamic pressure can be independently controlled by separately con- 
trolling the tunnel flow variables: total temperature, test-section mass 

flow, and the tunnel total pressure (ref. 3). The problem of closed-loop 
control of the tunnel total temperature, flow Mach number, and total pressure 
has been addressed and reported extensively in references 1 and 2. Utilizing 
the three tunnel control inputs, the liquid nitrogen valve area, the fan 
speed and the gaseous nitrogen bleed valve area, control strategies for 
precise closed- loop control of the tunnel variables have been derived and have 
been successfully demonstrated. 


1 Research Associate of Mechanical Engineering and Mechanics, Old Dominion 
University Research Foundation, P. 0. Box 6369, Norfolk, Virginia 23508. 
(on leave from National Aeronautical Laboratory, Bangalore-17, India). 


Since aerodynamic coefficients of test models are sought and analyzed 
as functions of the flow parameters: Reynolds number R g , flow Mach number 

M, and flow dynamic pressure q, a test program in a cryogenic tunnel 
obviously starts with the test-condition flow parameter set R , m, and q. 
Further, if the tests are not aeroelastic tests, the flow parameter set 
consists of R g and M only. The translation of the test flow parameter 
set R -M-q or R -M into the cryogenic tunnel variables P-T is the problem 
of test direction design. In this research effort, this problem of test 
direction design has been addressed, on the basis of minimum consumption of 
liquid nitrogen for the whole test program. 

In the foregoing analysis, a test direction parameter <p is defined. 

Its properties are used to derive ths order of test. Further, the problem 
of choosing the tunnel states P-T, such that minimum liquid nitrogen con- 
sumption occurs, is formulated and resolved. 

NOMENCLATURE 

a 
A 
c 
C 

P 

c 

v 
C 

m 
D 

E 

GN 2 
h 
J 
K 

LN 2 


velocity of sound in nitrogen (m/sec) 
area (m 2 or m 2 /m 2 ) 
chord (m) 

specific heat of nitrogen gas at constant pressure (kJ/kg-K) 

specific heat of nitrogen gas at constant volume (kJ/kg-K) 

specific heat of tunnel metal wall (kJ/kg-K) 

period of steady-state dwell (min) 

energy (kJ) 

gaseous nitrogen 

enthalpy (kJ/kg) 

joules 

gain constant or Kelvin 
liquid nitrogen 
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LNC 

• 

m 

mole 

M 

N 

P 

q 

R 

e 

R 

T 

t 

u 

U 

V 
W 

a 

6 

Y 


Newton 

* 

Z 

P 

u 


liquid nitrogen consumption 
mass flow rate (kg/sec) 

mole quantity CO. 028016 kg mass of nitrogen) 
test-section Mach number 
fan speed (rpm) 
pressure (atm) 

flow dynamic pressure (kg/m 2 ) 

Reynolds number 

universal gas constant (J/K/mole) 

gas temperature (K) 

time or time constant (sec) 

velocity (m/sec) 

internal energy (kJ) 

volume (m 3 ) 

mass (kg) 

cooling capacity of gaseous bleed (kJ/kg) 

cooling capacity of liquid charge (kJ/kg) = (121 + T) kJ/kg 

ratio of specific heats 

test direction parameter 

dot over a quantity refers to its time derivative 

1 kg mass accelerated at 1 m/sec 2 

aeroelastic test direction parameter 

compressibility factor (PV/RT) 

density (kg/m 3 ) 

viscosity (Newton-sec/m 2 ) 
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Subscripts 


a,b, . . .x 
1,2, . . .m, . . .n 
i 
s 

g 

F 

re 

m 

Q' 

low 

sat 

P 


integer indicies 

serial integer indicies 

integer index 

static 

gaseous 

fan 

Reynolds number 
motor 

consumption 
minimum usable value 
saturation 
pressure 


CRYOGENIC TUNNEL PROCESS 

A closed-circuit wind tunnel is a device intended for testing scaled 
models in fluid flow in order to determine their aerodynamic characteristics . 

The test-section fluid flow over the model is associated with flow parameters 
to which the aerodynamic characteristics of the model are related. For 
meaningful scaling of the model data to its full-scale value, it is necessary 
to maintain the flow similarity between the wind tunnel test-section flow 
and the full-scale flow. Amongst many flow parameters associated with flow 
similarity, Mach number and Reynolds number are the most important flow 
parameters for stationary models. In cryogenic wind tunnels, independent 
control of the flow parameters can be achieved because the tunnel conditions 
can be independently controlled (ref. 3). 

A closed-circuit cryogenic wind tunnel basically consists of a thermally 
autonomous pressure vessel designed as an aerodynamical ly contoured endless 
duct with a fan whose operation results in a range of flow velocities at the 
test section. The thermal autonomy of the cryogenic tunnel is realized by 
thermal insulation to the ambient either external to the metal walls or internal 
to the metal walls. The operation of the fan results in smooth motion of 
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the gas around the tunnel circuit accompanied by fan compression and wall 
friction heating. In order to maintain the tunnel gas temperature in the 
presence of this heating, liquified form of the test gas, liquid nitrogen, 
is sprayed into the tunnel. This controlled liquid nitrogen flow allows a 
cool down or regulation of the tunnel gas temperature. However, in this 
cooling process, the mass of the tunnel resident gas increases with 
accompanying increase of pressure. In order to maintain the tunnel pressure 
at a required value, it is necessary to remove a steady rate of warmer tunnel 
gas. This is usually achieved by controlled passive discharge of tunnel 
gas to the atmosphere. This complex interaction of fan-induced heating, 
liquid nitrogen induced mass increase and enthalpy decrease, and the gas 
discharge induced mass decrease and enthalpy decrease constitutes the 
cryogenic tunnel control problem. 

This thermodynamic cryogenic tunnel process has been extensively analyzed 
both graphically and analytically; an explicit, lumped parameter, nonlinear, 
multivariable model has been synthesized; a real-time, interactive simulator 
for the cryogenic tunnel has been developed, and finally the model and its 
real-time simulator have been validated by reconciling the steady-state, quasi- 
steady and transient responses of the 0.3-m Transonic Cryogenic Tunnel (TCT) 
to the simulator. This is presented in reference 1. 

The developed model has been used to analyze the closed- loop control of 
the cryogenic tunnel total temperature and total pressure. Further, two 
nonlinearly gain scheduled pID control laws have been developed with error 
magnitude dependent control logic to allow large set point imposition and 
smooth convergence to final set point. The temperature- loop control law 
has a feed-forward feature to allow rapid Mach number changes. These control 
laws have been successfully tested on the simulator and have been translated 
to a microprocessor-compatible flow chart. A microprocessor-based cryogenic 
tunnel controller has been custom built, and it allows fine tuning of the 
control process by allowing adjustments of up to 60 tunable coefficients of 
the control laws. This controller has been successfully commissioned on the 
0.3-m TCT to achieve temperature control to ±0.25 K and pressure to ±0.017 atm. 
The tunnel flow Mach number is directly controlled by the fan speed and can be 
run on a speed control loop or Mach number control loop. Details of this control 
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analysis, control law design, flow chart for microprocessor software and the 
successful operation of the controller are detailed in reference 2. 

Thus the problem of precise control of the cryogenic tunnel variables 
using the various control inputs has been solved. However, the problem 
of designing the set point set necessary to realize given flow has not been 
solved. In this report, the design of the tunnel set point set for desired 
aerodynamic test is addressed. 

The wind-tunnel user desires to perform tests on an aerodynamic model 
so as to obtain the characteristics of the model as a function of the flow 
parameters: Mach number M and, in the case of a cryogenic tunnel, Reynolds 

number R . The flow dynamic pressure becomes another relevant flow parameter 
if the tests involve aeroelastic behavior of the model. Otherwise, a 
minimal flow dynamic pressure q is usually sought. 

Translation of the flow parameters M, R , and occassionally q to the 

tunnel set points P, T, and N (or M) constitutes the tunnel test 

direction design problem. If all the flow parameters R , M, and q are 

a fixed set for a model, there exists a unique P, T, and N for each 

test. Then, the problem of test direction design is to determine in which 

order to perform a given set of tests. However, in most tunnel tests, the 

flow parameter test set R g and M are fixed with the constraint that the 

flow dynamic pressure q is to be minimum. Under these conditions, for 

every R e ~M combination there exists a range of P and T which can provide 

the necessary flow. However, by constraining the performance of the tunnel 

to minimum liquid nitrogen consumption or minimum q, an optimal combination 

of P and T can be determined for a given R and M. Under these 

e 

conditions, the problem of test direction design becomes one of determining 
the optimal P-T set and then ordering these for minimum energy consumption. 

FLOW PARAMETERS 
Introduction 

The fluid flow in the test section of a wind tunnel can be characterized 
by many dimensional and nondimensional flow parameters. The parameters of 
relevance to the present study are Mach number M, Reynolds number R g , and 
dynamic pressure q. 
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Mach Number 


The Mach number of a flow in the test section of a wind tunnel is 
related to the ratio of inertial forces to the elastic forces: 

^ _ u / inertial forces ^ 

a \ elastic forces 

The velocity of sound in the given a fluid is related to temperature and can 
be expressed as 


a = T (2) 

“mole s 

Reynolds Number 

The Reynolds number of a flow in the test section of a wind tunnel is the 
ratio of the inertial forces and viscous forces occurring in finite length 
of flow 


r - inertial forces _ puc ^ 

e viscous forces y 

The density of nitrogen gas in the test section of a nitrogen tunnel varies 
as a function of pressure and temperature (ref. 1). As a real gas this is 


P 


338.9 — (l + 250 — ^ 

rpi \ «p2 / 

s s 


kg/m 3 


(4) 


The velocity of the test gas as a function of gas temperature and flow 
Mach number is 


u 


20.38 JF 


M 


m/sec 


(5) 
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where, for nitrogen gas, a mole = 0.028 kg mass, y = 1.4, universal gas 

constant R = 8.31, and hence — = 20.38. 

mole 

The viscosity of the test gas, as a function of gas temperature, 
has been obtained by a simple function fit to Jacobson's data (ref. 5): 


U = 1.082 T s x 10~ 7 Newton-sec/m 2 


( 6 ) 


Hence the flow Reynolds number in a stream of gaseous nitrogen flow 
after appropriate approximation is 


- P 

R = £H£ = 65650 c M x 10 6 (7) 

e P T 1. 4 

s 


If the tunnel flow is isentropic, the following relationship between total 
and static quantities in the flow can be used with y = 1.4: 


P = P 
s 



; T s ■ T (l * M2 ) 


( 8 ) 


Then 


R 

e 


65650 — — c 

yl . 4 


M 

(l + 0.2 M 2 ) 2,1 


10 6 


(9) 
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Dimensionally, R^ has no dimension: 


puc 

u 


Newton sec 2 m 

■ ■ • ■ — — • 

3 m sec 

m J 


Newton - sec 



m 


( 10 ) 


For the 0.3-m TCT, which has a typical c = 0.1525 m. 


R = 10020 — 1 — 
e T i .4 


M 



0.2 


M 2 ) 


2. l 


10 6 


(ID 


In deriving the expression for Reynolds number of nitrogen flow, perfect 
gas behavior has been assumed in respect to density. However, since nitrogen 
gas behaves slightly differently, appropriate corrections for this can be 
made. Perhaps the best method for introducing this correction is to use 
the compressibility factor Z. Then 


puc 

e y 


P 

s 


ZRT 

s 


^ZYRT s 

y 


M c 


(12a) 


= 65650 P M - 1q6 

e .rr ti , i + / -A 2 * 1 

T (l + 0.2 M 2 ) 


(12b) 


However, for purposes of analysis, equation (9) has been used in subsequent 
treatment, and equation (12b) for final numerical computations as detailed 
in the Appendix (ref. 3) . 
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Dynamic Pressure 


Dynamic pressure is a flow parameter of structural loading significance 
and corresponds to stream kinetic energy transformed into potential energy 
as a structural load: 


q 



q = 


PM 2 y 

2 


10330 kg/m 2 


(13) 


Sensitivity Analysis 

The tunnel flow parameters R g and q are functions of the tunnel 
variables T and P with Mach number. The following sensitivity analysis 
provides an understanding of how the flow parameters vary with the tunnel 
variables : 


R 

e 


65650 — — c 
•pi . 4 


M 

(l + 0.2 M 2 ) 2,1 


The sensitivity of R g to P is 



The sensitivity of R g to T is 


(14) 


(15) 



The sensitivity of R g to M is 



(16) 


(17) 
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Similarly, the sensitivity of dynamic pressure q to tunnel variables can 
be determined. The sensitivity of q to P is 


5q 6P 
q - P 


(18) 


The sensitivity of q to M is 





(19) 


These sensitivity expressions can be used to determine the control accuracies 
and the sensor accuracies based on the desired flow parameter stability. 

In figure 1, the allowable error in the control of total temperature in 
a cryogenic tunnel for the given stability in Reynolds number of 0.5, 1, or 
1.5 percent is detailed. Typically the allowable error ranges from 0.35 K 
up to 3 K. 

In figure 2, the allowable uncertainty or error in the control of total 
pressure in a cryogenic tunnel for a given stability in Reynolds number 
of 0.5, 1, or 1.5 percent is illustrated.. Typically the allowable error 
ranges from 0.009 atm up to 0.1 atm depending on the tunnel pressure. 

It may be noted that the 0.3-m TCT closed-loop controller is able to hold 
temperature to ±0.25 K and pressure to ±0.017 atm, and hence is generally able 
to hold Reynolds number to about 0.6 to 0.7 percent. 

If the deviations in pressure, temperature, and Mach number occur 
concurrently, the uncertainty in Reynolds number is much more, and can be 
estimated as 



( 20 ) 
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By assuming the contributions of errors in P, T, and M to be weighted 
equally, for 1 percent uncertainity in R , the uncertainities allowed in 
tunnel conditions are about 0.4 percent in T, 0.6 in P and 0.6 percent 
in M. 

TEST DIRECTION PARAMETER 
Introduction 

A test direction parameter $ is now defined using the expression for 
Reynolds number in a cryogenic tunnel: 


R = 65650 c — 
e 


M 


1.4 


(l + 0.2 M 2 J 


2. 1 


10 6 


The test direction parameter <f> is defined as 


( 21 ) 



(l + 0.2 M 2 ) 


2. 1 


= 65650 c — — 10 6 

'j’l . 4 


( 22 ) 


The test direction parameter has unique properties which can be exploited in 
test direction design and economical operation of cryogenic wind tunnels. 

It can be observed that once the required flow parameters R g -M are known, 

4> is fixed and the tunnel temperature T and pressure P are uniquely 
related: 



(23) 
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where 


K = 65650 c x 10 6 
re 


Thus, knowledge of the test direction parameter uniquely describes choices 
of P and T as described by equation (23). For most cryogenic tunnel 
tests, the relationship of equation (23) provides the basis for optimal 
choice of the test direction. 

However, in aeroelastic tests involving fixed choice of R -M-q, the 
tests conditions P-T-M are fixed as in 


P 



1 

10330 


(24) 


and where T is defined as in equation (23) . 

Properties of Test Direction Parameter 

The test direction parameter <f > uniquely fixes the ratio P/T 1,4 . 

Using the 0.3-m TCT as an illustrative example, a plot of R g vs. M for 
loci of constant <j> is shown in figure 3. The tunnel flow Mach number 
varies typically from 0.4 to 1, and the flow Reynolds number from about 
0.5 x 10 6 to 72 x 10 6 . Eleven loci of <f> varying from 5 x 10 6 to 105 x 10 6 
in steps of 10 x io 6 are shown in figure 3. 

Consider a typical constant $ locus A 2 E 2 corresponding to <J> = 105 x 10 6 
This locus represents a range of Reynolds number from about 39 x 10 6 at 
M = 0.4 to 72 x 10 6 at M = 1.0. The locus A 2 E 2 corresponds to fixed P/T 1 - 4 . 

The choice of tunnel total pressure P and temperature T of a cryogenic 
tunnel is limited to its operating envelope. Figure 4 illustrates the tunnel 
ideal operating envelope 1 <_ P 6 atm and T^^ <_ T <_ 350 K. The pressure 
envelope limit is imposed by the maximum safe pressure based on appropriate 
pressure vessel code and by the fact that for passive gas bleed the tunnel 
must function at or above atmospheric pressure. The low temperature limit is 
imposed by the need to keep the tunnel flow in gaseous phase. The upper 
limit is dictated by structural considerations. 
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Figure 4 provides the loci of constant test direction parameter <J> in 

the plane P vs. T. Eleven loci of <j> = 5 x io 6 to 105 x 10 6 in steps of 

10 x 10 6 cover the tunnel envelope. Each locus corresponds to the choice of 

P-T combination for a given R g -M and hence <j>. Consider a typical locus 

DE corresponding to <j> = 15 x io 6 . A choice of P from 1 to 6 atm and T 

from 115 to 350 K can be made along line DE. Further, <p corresponds to a 

range of R g from 5.6 x 10 6 to 10.2 x io 6 and of M from 0.4 to 1.0. 

Thus figure 4 provides an indication of the flexibility involved in translating 

R -M to P-T. 
e 

Figure 4 further illustrates the limitations imposed by the tunnel fan 
power and fan speed derived using the following identities from reference 2: 


K P^T M 3 

Power limit = 1900 kW = — — ('25') 

(l + 0.2 M 2 ) 3 

-0.035 

Fan speed limit = 5600 rpm = 597 M (1 - 0.3 M) P (26) 


The power and speed boundaries limit the choice of P and T for a given 
M and cf>. 

It is evident from figure 4 that, at a given <p, the tunnel total 

pressure increases along the locus DE and hence the flow dynamic pressure 

q also increases along the locus DE. Although the minimum ideal tunnel 

total pressure boundary is 1 atm and the minimum ideal total temperature 

boundary is T the minimum usable total pressure in passive 

bleed mode is about 1.2 atm and the minimum usable total temperature is 

a T low' w hich exceeds ^ sat by a safe margin. The choice of lowest 

usable tunnel pressure for a given <p occurs along the locus AB^ B 

of figure 4, and either along the T locus BC or along the safe 

s a l 

temperature locus T low corresponding to BjCj. Hence the ideal locus 
ABC or the safe usable locus A3jCi turns out to be the minimum dynamic pres- 
sure locus for any set of test <t>. As will be demonstrated later, the locus 
ABC/AB^Ci a lso turns out to be the minimum liquid nitrogen consumption 
locus for a given set of test <f>. Thus figures 3 and 4 provide an insight 
into the properties of test direction parameter <}> in respect to R g -M and 
P-T as well as to q. 
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MINIMUM ENERGY TEST DIRECTION DESIGN 


Introduction 

Consider a model test program in a cryogenic wind tunnel wherein aero- 
dynamic data on the model is desired at a number of flow parameter conditions 
of R^-M combinations. Depending upon the number of angles of attitude, a 
minimum dwell time of steady flow is associated with each flow condition. 

Let the flow parameter condition set be expressed as (Rg-M-D)^, (r^-M-d)^, 

. ..(r^,M,d)^ where subscripts a,b,...x correspond to arbitrarily ordered 
n conditions of test, and D corresponds to dwell time. 

The problem of test direction design in a cryogenic wind tunnel is, 
firstly, the translation of the flow parameter test condition ^r^-M-d)^ 
to the tunnel variable condition (p-T-M-d)^ and, secondly, ordering the 
integer set i = a,b,...x to an execution order i = l,2,3,...n such that overall 
performance is optimal either in terms of minimum liquid nitrogen or least 
dynamic pressure or both. 

The problem of determining the minimum liquid nitrogen consumption 
strategy in realizing a given test program in R^-M is now considered in two 
cases. In the first case, the liquid nitrogen consumption in cooling the 
tunnel shell is ignored and the problem of choosing the best P-T combination 
for least liquid nitrogen consumption is addressed. In the second case, the 
total liquid nitrogen consumption in cooling the tunnel, in buildup of 
tunnel pressure and in the best choice of P-T for given R^-M such that 
least consumption occurs is detailed. 

Case 1: Minimum Liquid Nitrogen Consumption Strategy, Considering Fan Heat 

Only, for a Given Test Program in R^-M 

The fan energy released into the tunnel gas stream when the tunnel 
conditions are pressure P, temperature T and Mach number M is 

K_ P JT M 3 

E = — ! - - : 60 D (27) 

(l+0.2 M 2 ) 
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The liquid nitrogen consumed in maintaining the tunnel temperature at 
T is (ref. 1 ), 


LCN = -=■ = 


E 60 D K p P yfT 


M 3 


6 (l + 0.2 M 2 ) 3 (l21 + t) 


kg 


(28) 


In a given test program, with the desired set of R -M known, the test 

P e 

direction parameter <P and hence — is fixed. As described by equation 


(23): 


1 

1.4 


T = [ P 


' 1. 4 


P = q> T 1 

K 

re 


where 


K = 65650 c x 10 6 
re 

Substituting this expression in equation (28) , the liquid nitrogen consumption 
LCN over a period of 60 D seconds of steady flow can be expressed as a 
function of temperature only: 


LCN 


T l . 9 

K Q 121 + T 


(29) 


where 


K _ 60 D <fr M 3 

Q (1 + 0.2 m 2 ) 3 k 

' re 

The liquid nitrogen consumption can be minimized as a function of T: 


a (LCN) 
3(T) 


0.9 

T 


1.9 (121 + T) 

(121 + t) 2 



(30) 


16 


1 


! 




\ 


I 

! 


-I » 

i 


Hence extrema occur at T = -255 K or T = 0. Obviously these extrema do not 
occur in the physically realizable range of T < T < 350 K. The lowest 
value of LCN occurs at the lowest physically realizable value of T = T gat 
as evident from inspection of equation (29) . 

However, in operating a cryogenic tunnel, the stream local static 
temperature over the model is far less than the stream total temperature T. 
Hence a safe margin is necessary in choice of T^ ow such that T^ qw = T sat 
+ margin. The magnitude of margin is a function of flow Mach number. In 
reference 6, determination of safe margin based on theoretical and experimental 
analysis is detailed. For purposes of this report, a safe maximum margin of 
12 K has been used for all computations. The saturation vapor phase boundary 
of nitrogen can be expressed, on the basis of simple function fit for data 
from reference 5, as 

T = 50 + 27.34 P °- 296 kelvin ( 31 ) 

sat ' 

The locus of vapor phase boundary T is shown in figure 4 as locus BC. 

S 3. u 

The liquid consumption can also be expressed as a function of pressure P 
alone as 


LCN = K 

P 


1 . 36 
P 


0.71 

K x + P 


This function can be minimized for P: 



(32) 


(33) 


It can be observed that again the minima occur at either P = 0 or 
P = -fifKj). Both these values of P are physically unrealizable since they 
fall out of the tunnel range P^ QW £ P <_ 6 atm. Further, it can be noted that 
lowest LCN occurs at P^ ow i- n the physical world as can be observed from 
equation (32). 
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Hence, given a set of R g -M combinations to be realized, the best 
operating P-T combination occurs either along the locus P^^-T or 
Tl 0w -P locus. The least liquid nitrogen locus ABjCi, wherein ABj is the 
P iow -T segment and B^C^ is the T^^-P segment, is shown in figure 4. 

Obviously, given the desired R g -M combination, P^ QW is assumed and T 
determined as per equation (23) . The point A corresponds to the lowest 
<J and point C corresponds to highest p. In the 0.3-m TCT, segment 
AB corresponds to p from 1 x io 6 to 27.5 x 10 6 . Segment BC corresponds 
to <p from 27.5 x 10 6 to 105 x 10 6 . However, the usable locus B^Ci 
corresponds to p from 24 x io 6 to 87 x io 6 . 

The problem of ordering the flow parameter test set (R e -M-D) a , 

d)^, ... (R e -M-D) x is now considered. Each operating flow parameter 

point is associated with a finite p, and hence the test program set has 

a test direction parameter set of <f> ,<J>, ,...<(> . The numerical values of 

p ,p, ,...p are determined, and these are reordered such that < Pz <••• 

T a r b’ T x ~ — 

<j>^. The test direction providing least liquid consumption is <f>i ,p2> • • • ( t , n - 

Consider the first test point 4> x , corresponding to (R g -M-D ) 1 belonging 
to the set 4) ,...<{> . This test direction parameter <{>1 is now translated into 

3. X 

the tunnel operating condition Pi,Ti,Mj,Di . Since the value of i is the 
lowest, as shown in figure 4, the least energy operating point is likely to 
be along the locus AB^. The choice of pressure is Pj = P^ QW = 1-2 atm, and the 
complementing temperature Tj is determined from equation (23), in the 
range T^ qw <_ Tj <_ 350 K. After performing tests at <f>i for period Dj the 
tunnel is moved to Pz- Note that Pz i • 

The new tunnel pressure ?z and temperature T 2 are estimated from 
equation (23), at the low tunnel pressure of P 2 = p i = P low > Obviously, 
temperature T 2 <_ T x in the range T 1qw < T 2 < 350 K. In reaching T 2 , cooling 
of the tunnel shell from T x to T 2 is involved. This process is carried 
out for ij> 3 , . . . up to a <t> i wherein T i < T 1qw> When this stage is reached 
the new tunnel conditions are estimated on the basis of T^ = T^ qw , and P^ is 
estimated from equation (23). Subsequently, + 1 to are trans lated 

to P^ + i - T^. + 1 , etc. and are likely to occur along locus B 1 C 1 of figure 4. 

In this process of translation of [<j>] to [P-T], it is obvious that, 
since 44 < Pz <_. . - <t> n the tunnel conditions also bear the relation ?i P 2 i p 3--* 
> P n and Ti <_ T 2 <_. . . along locus ABj. However, the tunnel temperature bears 
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the relation T > T . >... along locus B 1 C 1 , since the vapor phase boundary 
n — n-1 — 

has positive slope in the P-T plane. Thus the execution order involves no 
reversals in tunnel conditions and proceeds monotonically from near ambient 
P-T to tunnel extremes, thereby assuring minimal consumption of liquid 
nitrogen. 

Case 2: Minimum Liquid Nitrogen Consumption Strategy, Considering Fan Com- 

pression, Metal Shell Cooling, and Mass Buildup, for a Given Test Program in 

R -M 
e 

Consider a cryogenic wind tunnel inactive at ambient conditions. To use 
this tunnel to realize a tunnel flow of R g -M for a period of 60 D seconds 
and then terminate the test, a finite amount of liquid nitrogen consumption is 
involved. The problem of minimizing this total liquid nitrogen consumption 
as a function of temperature is now considered. 

The energy required to cool the tunnel from ambient of say 300 K to a 
temperature T is dominantly controlled by metal enthalpy between the 2 
temperatures. This has been experimentally verified for the 0.3-m TCT 
(ref. 1). To an acceptable degree of accuracy, the liquid nitrogen required 
to cool the tunnel from 300 K to T, at any pressure is 


w r + w c 

8-51 (300 - T) kg ( 34 ) 

8 

It is assumed that tunnel cooling occurs at lowest pressure of 1.2 atm and a 

low convection Mach number of M = 0.3 or less. In estimating the^mass of gas 

W , a pressure of P can be used at an average temperature of 5 * 

8 , 

The liquid nitrogen consumption in building the tunnel pressure can be 

estimated from the density expression at constant temperature as 


P - 1 

341 V kg 

where the pressure buildup is from 1 atm to P. 


(35) 
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The total liquid nitrogen consumption in performing a cryogenic tunnel 

test at R -M for 60 D seconds is 
e 


60 


LCN = 


D K p PVT M 3 


W C + W C 
g v t m 


(l + 0.2 M 2 ) 3 (121 + T) (121 + T) 


(300 - T) + (P - 1) 


(36) 


Since the test direction parameter $ is known from equation (23) as 



Hence 


LCN 


1.9 


(121 + T) 


W C 
t m 


W C 
_g_v 


(121 + T) 


(300 - T) + 


341 V 



(37) 


The liquid nitrogen consumption consists of three terms. As temperature T 

moves down from 300 K, the second term grows linearly. At its maximum, 

when T = T, the second term has a maximum of about 1800 kg. The third term 
low 

grows linearly with pressure, and, at the maximum tunnel pressure of 6 atm, the 
liquid nitrogen involved corresponds to about 350 kg. The first term is 
power and dwell period D dependent. 

It can be noted that cooling the tunnel involves nearly 1800 kg of 
liquid which can be eliminated if the tunnel run can be performed at high 
temperatures. The minimization of equation (36) analytically with respect to 
T would yield this solution. In the present analysis, a numerical estimation 
of LCN has been carried out as a function of temperature and the results are 
presented in figures 5 to 14. These correspond to R g of 10 x 10 6 to 
50 x 10 6 in steps of 10 x io 6 at M = 0.6 and M = 0.9. The program generates 
combinations of P and T necessary to realize given R^-M, while conforming 
to physical realizability. 


20 


Figure 5 is a plot of LNC as a function of T in the pressure range 
of 1.2 < P <_ 6 for various dwell times ranging from 2 to 22 min. When the 
dwell period is 2 min, the test program for realizing R g = 10 x 10 6 at 
M = 0.6 is most economical only at 300 K and not 105 K. However, as the 
dwell time increases, the loci of LCN which have a negative slope at D = 2 
gradually turn positive. At a 22-min dwell time, the locus has a positive 
slope, indicating that operation of the tunnel at 105 K is most economical. 
Obviously the metal enthalpy term of 1800 kg is more than saved at high 
dwell times and hence makes cool down worthwhile. 

In figure 6, R g = 20 * 10 6 at M = 0.6, the range of usable temperature 
is reduced. The issue here is whether to cool the tunnel to 200 K or down 
to 85 K. Again the loci of LCN provide the answer in that for low dwells 
full cool down is not necessary, whereas, for large dwell periods in 

excess of 10 min, it is necessary to cool the tunnel down to save liquid 

nitrogen. 

The progressive reduction of temperature range for R g = 30 x 10 6 to 

50 x 10 6 at M = 0.6 can be noted in figures 7 to 9. In figure 9, the choice 

of temperature is only a few degrees. 

Similar results are illustrated in figures 10 to 14 for R g = 10 x 10 6 
to 50 x 10 6 at M = 0.9. Again, similar conclusions regarding dwell time can 
be observed. In summary, as long as the cryogenic tunnel experimental dwell 
time is of considerable length, involving a reduction in liquid consumption 
by cool down in excess of 1800 kg, the least energy consumption locus is 
the minimum temperature/minimum pressure locus. However, for very short 
experiments there exists another choice of operation at highest feasible 
temperature, with the associated higher q flows. Major cryogenic tunnels 
are not likely to run for such short periods, and so the second optimal 
solution is meaningless. 


In the previous sections, the problem of test direction design in a 
cryogenic wind tunnel was defined as one, firstly, of ordering the desired 


LOCUS OF OPTIMAL TEST DIRECTION 


flow parameter conditions 
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and, secondly, of translating these operating flow parameters into tunnel 
conditions (P-T-M-D) i , . . . (P,T,M,D) n . In doing so the constraint of 
minimizing liquid nitrogen consumption was fulfilled. In addition, the 
chosen flow dynamic pressure was the minimum ever achievable. 

The best tunnel operating locus in the tunnel pressure-temperature 

envelope turned out to be along the minimum usable tunnel pressure line at 

varying temperature in the range T^ qw <_ T <_ 350 K and secondly along the 

minimum usable temperature line at varying pressure in the range P^ ow £ P 

<_ 6 atm. The locus AB^C^ of figure 15 illustrates the best operating locus. 

Surprisingly, the inference that can be drawn by this study is that a 

cryogenic tunnel need not operate at all at high pressures when at any 

temperature other than T^ ow * The locus ABjCi is capable of providing the 

total performance capability of the cyrogenic tunnel in R -M. Need for 

e 

working at any other point in the P-T plane of figure 4 or 15 arises only 
if R -M-q combinations are desired such that q is a dynamic pressure higher 
than the minimum along line AB^Ci. 

The best test direction locus ABC is transformed into the energy state 

plane U vs. W and is shown in figure 16. The energy state diagram shows the 
§ 

internal energy-mass plane with loci of constant pressure, constant temperature, 
and constant test direction parameter <f>. The minimum liquid nitrogen consump- 
tion test direction locus ABiCi moves, firstly, along the minimum tunnel pressure 
locus ABi and then along the minimum temperature locus B^Ci. 

The test direction design problem can be exemplified by points A 2 , B 2 , C 2 , 
and D 2 of figure 3. Consider four Reynolds number-Mach number flows of 
R g = 39 x 10 6 - M = 0.4 at point A 2 , R g = 42 x 10 6 - M = 0.58 at point b 2 , 
r = 48 x 10 6 - M = 0.57 at point C 2 , and R = 51 x io 6 - M = 1.0 at point d 2 . 
The normal tendency for a tunnel user would be to perform tests in order 
A 2 B 2 C 2 D 2 . However the estimation of <)> for each case shows that is 
105 x 10 6 , 85 x io 6 , 95 x 10 6 , and 75 x 10 6 , respectively, for the points 
A 2 , B 2 , C 2 , and D 2 . The preferred order of execution as analyzed previously 
turns out to be D 2 B 2 C 2 A 2 in ascending order of <P and would involve tests 
during a pressure buildup at Tj . 
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The problem of test direction design in a test program wherein the 
desired conditions are (R e -M-D-q) a> (R e ~M-D-q)^, . . . (R e ~M-D-q) n can now be 
considered. This type of test would probably be an aeroelastic test. 

Since all the three flow parameters are given, there is no minimum liquid 
consumption problem. However, the problem of ordering does exist theoretically. 
Given (r -M-D-q), the choice of P-T is unique, since q = — 10330 



1 

10330 


(38a) 


Then 


T = 



0.2 M 2 ) 


2.1-1 


1 . 4 


65650 c M 


(38b) 


The problem of ordering could be conceived through another test parameter \p. 




+ 


2N 2 * 1 
0.2 M 2 ) 

~q~ 


a MZ i 0 6 

■p 1 . 4 


(39) 


This expression can be easily derived from equations (9) and (24). 

The dynamic pressure ranges from 1,168 kg/m 2 at P = 1 atm and M = 0.4 to 
43,400 kg/m 2 at P = 6 atm and M = 1.0. The test direction for aeroelastic 
tests can be ordered in <|». However, since ip is directly proportional to 
•p" 1 • 4 ^ the aeroelastic tests can be planned monotonically without having to 
conceive of or evaluate f. Hence, for aeroelastic tests, equations (38a) and 
(38b) are used to evaluate P and T, and then the tests are directly ordered 
in T. Thus (R e -M-D-q) a ,... (R e -M-D-q) x is transformed to (p-T-M-d) s , .. . 
(p-T-M-d) . The temperature values T & , — T^ are then used to reorder them 
such that Ti > T 2 >...T where l,2,...n are reordered set a,b,...x. 
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TYPICAL CRYOGENIC TUNNEL TEST DIRECTION DESIGN 


Introduction 

A computer program in Fortran which provides the typical test direction 
design for a given test program consisting of desired Reynolds number-Mach 
number flow conditions is detailed in the Appendix. The program is a simple 
mechanization of the test direction design procedure analyzed in previous 
sections. The program uses fixed tunnel parameters specific to the 0.3-m 
TCT viz volume of 14.1 m 3 , metal structure weight of 3200 kg, fan power 
constant, and the cooling capacity of liquid nitrogen. 

The program initially reads in the number of test points n and 
associated flow parameter set R^-M-D. The test direction parameter <f> is 
determined for each of the flow parameters, and the test flow parameter 
set is reordered such that <f>i £ 4>2 _i • • • • Secondly, the program translates 

each flow parameter condition R^-M-D into tunnel test condition P-T-M-D 
such that minimum liquid nitrogen is consumed in the process. The estimate 
of liquid nitrogen required for fan power cancellation and for transition of 
the tunnel from the previous condition to the present condition is estimated. 

It is assumed that the metal shell cooling occurs at low Mach number and low 
pressure, thereby assuring minimum wastage of liquid. The computer program 
incorporates real gas correction in the form of compressibility effects. 

Example 

A specific test program realizing 38 flow parameter combinations ranging 

from R = 0.9 x 10 6 at M = 0.1 to R = 50 x 10 6 at M = 0.8 is now considered, 
e e 

The table in figure 17 describes the 38 flow parameter conditions with the 
estimated dwell time. 

The problem is to decide the order in which these flow parameter conditions 
are to be realized in the cryogenic tunnel and to translate the flow parameters 
to tunnel conditions P-T which provide minimum consumption of liquid nitrogen. 

The program is used to evaluate the test direction parameter <f> in the 
order <p i _< <f> 2 - - 4>38 • The result of this process is shown in figure 18. The 

tests in the tunnel are to be performed in the same order. Note that the 
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lowest R = 0.9 x 10 6 at M = 0.1 is listed as the twelfth test point and the 
6 

highest R g = 50 x 10 6 at M = 0.84 as the thirty-first test point. 

The Fortran program has been used to translate the flow parameter test 
conditions of figure 18 to the tunnel operating conditions P and T. The 
program further estimates the liquid nitrogen consumed to cancel fan power 
and liquid nitrogen consumed for tunnel condition transition. The tolerances 
on the tunnel pressure and temperature control to assure half percent 
uncertainty in the flow Reynolds number are also evaluated. The table of 
figure 19 thus provides the solution to the problem of test direction design 
for conditions tabulated in figure 17. 

It may be noted that the first test point starts at the lowest tunnel 
pressure of about 1.2 atm and the highest tunnel temperature of 235 K. The 
subsequent operating points occur in the cool down direction. When the 
test <p reaches 24 x 10 6 , corresponding to R g = 15 x 10 6 at M = 0.8 
the temperature limit is encountered. Subsequent test points occur at the 
low temperature limit of T ^ + margin, but at increasing pressures. It is 
obvious that the tunnel conditions move monotonically through the test range, 
assuring minimal transition in liquid nitrogen consumption. 

The test program can be run in the reverse order also without any extra 
loss of liquid nitrogen. The program finally prints the total liquid nitrogen 
consumption estimate to perform the desired test. 

CONCLUDING REMARKS 


The cryogenic wind tunnel is a device intended for generating very high 
Reynolds number flows past aerodynamic models. This increase of flow Reynolds 
number is obtained by cooling the test gas, nitrogen, by use of liquid nitrogen 
sprayed into the tunnel. 

The cost of using a cryogenic tunnel is dominated by the cost of liquid 
nitrogen. Typically the cost of energy in the form of liquid nitrogen is 
nearly 80 to 100 times the cost of equivalent electrical energy. Thus the 
liquid nitrogen consumption in a cryogenic tunnel test program should be 

minimal . 
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An important requirement of the cryogenic tunnel use is that the user 
should conceive of the wind tunnel test program on the basis of desired flow 
R g } flow M, and dwell period. The dwell period depends upon the number 
of attitudes of the model and data scan time. In the contemporary, ambient 
temperature tunnel usage practice, the tunnel users usually conceive of the 
test on the basis of Mach number and tunnel blowing pressure. 

In this report, analysis of techniques for ordering a test program 
and translating the tunnel flow parameters into the tunnel flow conditions 
P and T has been carried out. A test direction parameter <j> has 
been defined to provide the basis for ordering the test program. Secondly, 
a simple extremum determination analysis has been carried out to determine 
the best tunnel conditions P-T involving least energy consumption. This 
problem has been shown to be valid only when the flow dynamic pressure q is 
unconstrained. Further, the test liquid nitrogen consumption locus has been 
shown to be the least dynamic pressure locus also. 

In the case of aeroelastic tests involving R e ~M-q, there exists a 
single P-T-M set for each flow condition. The test order here is based on 
tunnel temperature T being chosen monotonically. 

In the procedure used for translating the flow parameters to tunnel 
variables, a single real gas correction for compressibility has been made. 
Better accuracy can be obtained in translating R -M to P-T by incorpo- 
rating viscosity effects. However, this does not affect the test direction. 

In conclusion, the problem of test direction design in cryogenic wind 
tunnels has been defined, analyzed and basically resolved in this report. 
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APPENDIX 


PROGRAM LISTING FOR TEST DIRECTION DESIGN 


PROGRAM MAIN (INPUT* OUT PUT,TAPE1, TAP E5- INPUT, TAPE6- OUT PUT) 
DIMENSION M(100),MC!100),U!100)>RE(100),D!100),R(100),UR<100), 
IP ! 100 ) , T ( 100 ) » DL'M { 100 ) » Dt P ! 100 ) » OET ( 100 ) 

CCC PROGRAM VARIABLES 
C INPUT FLOW PARAMETERS 
C RE •REYNOLDS NUMBER * 10E-6 
C M ■ MACH NUMBER 

C U • TEST DIRECTION PARAMETER * 10E-6 
C D • DWE LL> M INUTE S 
C EXECUTION ORDER FLOW PARAMETERS 
C P • REYNOLDS NUMBER * 10E-6 
C MC • MACH NUMBER 

C UP - TEST DIRECTION PARAMETER * 10E-6 
C OLM • DWELL, MINUTES 
REAL M,MC 
N»3e 

CCC ARC- MARGIN IN TEMPERATURE TO ONSET OF CONDENSATION 
ARG-12 

CCC N • NUM8ER FLOW PARAMETER TEST POINTS IN PROGRAM 
G-2.1 

CCC READ INPUT DATA- FLOW PARAMETER SETS 
READ 4, (RE(I),M(I),D(I),I-1,N> 

4 FORMAT! 3F 10. 2) 

PRINT 16 

16 FORMAT! 1H1, / / 1 2X+ R E* 1 3 X *M* 1 3X *D* ) 

PRINT 17, ( PE ( I ),M ( I ),D( I), 1*1, N) 

17 FORMAT ( 3F14 . 2 ) 

CCC DESIGN TEST ORDER BY MONOTONIC ORDERING OF TEST DIRECTION 
CC PARAMETER U- R /M * ( ( 1 + 0 . 2*M*M ) ♦* 2 . 1 ) 

DO 5 J«1»N 

5 U(J)-RE!J)*((1+0.2*M(J)*MCJ) )**G ) /M ( J ) 

DO 6 J ■ 1 , N 

UU--0.001 
DO 10 1*1, M 
X-U! I)-UU 

IF(X.GT.O) GO TO 15 
GO TO 10 

15 

L « I 

10 CONTINUE 
NN-N-J+1 
R ( NN ) *RE ( L ) 
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MC (NN)-M(L) 

UP (NN )-U(L) 

DHM(NN)*0(L) 

U( l ) *-C.C01 
6 CONTINUE 
PPINT 19 

19 F0PMAT(1H1»//6X*DIRECTIQN*8X*PE*13X*M*) 

PRINT 20 1 (UR (I ),R ( I ) # MC (1 ) , I - 1 # N ) 

2C FOPMAT(3F14.2) 

PRINT 22 

2? F0PMAT(1H1W/19X*MARGIN TO ONSET OF CONDENSATION^*) 

PPINT 2 3 » AR G 
23 FORMAT! 31X, F4.1) 

PPINT 21 

21 FGPMA T ( //5X*PE*9X*M*7X*P*8X*T*6X*DWELL*4X*FANLN*5X*TLN* 
14X*DE°*3X*DET* ) 

CCC TRANSLATE RE-M IN TO TUNNEL CONDITIONS P-T 
SI'M-0 
PPP-1 
TTT* 3 00 
DO 100 1*1, N 

CCC START AT PRESSURE OF PLOW = 1.2 FOR 0.3-M TCT 
P ( I ) * 1 • 2 

25 CONTINUE 

TT»< P < I )*10020*MC < I ) /P < I ) / ( { 1+0.2*MC ( I) *NC ( I > )**G) >**( 1/1.4 ) 

CCC ADJUST FOR COMPRESSIBILITY FACTOR 

26 CONTINUE 

DC*1. 37-e.773E-2+TT+4.703E-4*TT*TT-l. 336E-6* ( TT**3 )* 
11.462E-9*(TT**4) 

EE«5,521— 1.986E-1*TT+7.817E-4*TT*TT-1.253£-6*(TT**3)+ 
15.333E-10*(TT**4) 

ZZ*1-EXP<D0)*P!I>-EX?(EE)*(P< I)**2> 

IF (ZZ.LE.O. 5 JGO TO 31 
GO TO 32 

31 ZZ-0.5 

32 RP-10020 + P! I ) / (TT**1.4 )*MC < I) /( (1*0.2*MC( I )*MC! T ))**2.1 ) /SORT <ZZ) 
CCC ITERATE FOR COMPRESSIBILITY ADJUSTMENTS TO REYNOLDS N U M 3 E P 

DR* ABS ! ( R ( I ) -P P ) / R ( I) ) 

IF(DR.LE. 0.002 JGO TO 27 
IF(R( I) .GT.RR)GO TO 28 
IF t R ! I ) • LT.RR ) GO TO 29 
28 TT *TT-0 .01 
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GC TO 26 

29 TT-TT+O.C1 
GP TO 26 

27 CONTINUE 

CCC CHECK FOR PHYSICAL REALISABILITY OF TEMPERATURE 
CCC CHECK FOP ONSET OF CONDENSATION 

CCC VAPOUR PHASE BOUNDARY IS 50 + 2 7 . 3 A * ( P *70 . 296 ) 

X«TT-50-27.3A*(P< I ) **0 . 296 ) -ARG 
IF(X.LT.O) GC TO 30 
GO TO *0 

30 P( I ) «P( I )+0.01 
GO TO 25 

AO T ( I ) -TT 

R1*R ( I ) 

P2«P< I ) * 1 . 005 

P1-P1*(T(I)**1.A)*((1+0.2*MC ( I ) *MC ( I ) )**G) /10020/MC ( I ) 
P2»R2*(T(I)**1.A)*((1+0.2*MC(1)*MC(I) )**G ) /10020/MC ( I ) 

DEP( I )■( P1-P2)*1A .696 

T1“(°1*10020*MC(I ) / R 1 / < ( 1 + 0.2*MC ( I ) *MC ( I ) )**G ) ) ** ( 1 / 1 . A ) 
T2«(P1*1C020*MC(I)/R2/((1 + 0.2*MC ( I)*MC( I) ) **G ) ) ** (1 / 1 . A > 

DET( D-T2-T1 

CCC CHECK FOR PHYSICAL R E AL I S AB I L I T Y OF PRESSURE 1.2<P<6 
Y»6-P ( I ) 

IF (Y.LT.O) GO TO A5 
GC TO A 7 
A 5 CONTINUE 
C-0 TO 100 

CCC ESTIMATE LN2 REQUIRED FOR FAN COMPRESSION HEAT CANCELLATION 
A7 H-57.2*P(I)*SQPT(T(I))*(MC(I)**3)/(121 + ni)) 
H-H/((1+0.2*MC(I)*MC(I))**3) 

S * H*DUM ( I ) *6 0 
X X »TTT-T ( I ) 

IF(XX.LT.0)G0 TO A3 
GO TO 44 
A3 TTT-T(I) 

AA CONTINUE 

CCC ESTIMATE LN2 REQUIRED FOR TUNNEL SHELL COOLING £ PRESSURISATION 
S1*(TTT-T(I) )*9+( P( I ) - P P P ) *A779 / I ( I ) 

SUM-SUM+S+S1 
PPP-P ( I) 

TTT-T(I) 

CCC PRINT EXECUTION ORDER * TUNNEL STATES, DWELL TIME, ESTIMATE D LN2 CONSUMPTION 
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C CONTROL ERROR TOLFRANCES ETC 

CCC DET ■ TOLERANCE 1^ TEMPERATURE# K FOR 1/2 PERCENT STABILITY IN RE 
CCC 0 E P ■ TOLERANCE IN PRES$URE#PSI FOR 1/2 PERCENT RE STABILITY 
C TRAN LN2 * LN 2 USED FOR SHELL CUOL DOWN £ P R E S SUR I S AT ION 
A8 PRINT 50> (R ( I ) #MC ( I ) # P { I ) # T ( I ) , DUM( I ) #S# SI# DEP (I )# DFT ( I ) ) 

50 F0R M AT(4F9.2#3F9. 1,2F6.2) 

100 CONTINUE 
PRINT 59 

59 FORMAT ( //2Y*T0T4L L N2 USED*) 

PRINT 60# SUM 

60 FORMAT (F12»3) 

STOP 

END 
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Figure 1. 


Allowable error in control of tunnel total temperature for a 
given stability in Reynolds number in a cryogenic wind tunnel. 
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Figure 2. Allowable error in control of tunnel total pressure for a 

given stability in Reynolds number in a cryogenic wind tunnel 
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Test direction parameter =_!< l + 0.2 M 2 ) 1 1 = 69650 — ~ 10 6 
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Figure 3. Test direction design in cryogenic wind tunnels — test direction 
parameter 4> as a function of R g and M. 
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Test direction design in cryogenic wind tunnels— locus of least 
liquid nitrogen consumption and least q in P-T plane. 
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Figure 5. Liquid nitrogen consumption in a cryogenic tunnel test at 

R = 10 x 10°, M = 0.6, and D = 2 to 22 in 4-min steps, 
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Figure 6. Liquid nitrogen consumption in a cryogenic tunnel test at 
R = 20 x 10 6 , M = 0.6, and D = 2 to 22 in 4-min steps. 
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Figure 7. Liquid nitrogen consumption in a cryogenic tunnel test at 
R g = 30 x io 6 , M = 0.6, and D = 2 to 22 in 4-min steps. 
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Figure 8. Liquid nitrogen consumption in a cryogenic tunnel test at 
r = 40 x 10®, M = 0.6, and D = 2 to 22 in 4-min steps, 
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Figure 9. Liquid nitrogen consumption in a cryogenic tunnel test at 
R e = 50 x 10°, M = 0.6, and D = 2 to 22 in 4-min steps. 
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Figure 10. Liquid nitrogen consumption in a cryogenic tunnel test at 

R = 10 x 10*, M = 0.9, and D = 2 to 22 in 4-min steps, 
e 


41 


R * 20 10 6 M • 0. 9 0 • 2. 22. d mts 
s 


Fan power 

LNC 60 0 

( 121 ♦ T ) 


( 300 -T )■ 


W. C + W. C 
l m <3 v 

( 121 ♦ T ) 


341 V 

( ?- 1 ) — 



Figure 11. 


Liquid nitrogen consumption in a cryogenic tunnel test at 

R = 20 x 10 , M = 0.9, and D = 2 to 22 in 4-min steps, 
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Figure 12. Liquid nitrogen consumption in a cryogenic tunnel test at 
R ^ 30 x 10®, M = 0.9, and D = 2 to 22 in 4-min steps, 
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LNC is the LNj used to cool tunnel 

from 300K to T, to develop pressure 
from 1 atm to P and to provide flow 
at R g & M for 60 0 seconds. 
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Figure 13. Liquid nitrogen consumption in a cryogenic tunnel test at 
R g = 40 x 10°, M = 0.9, and D = 2 to 22 in 4-min steps. 
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Fieure 14. Liquid nitrogen consumption in a cryogenic tunnel test at 
R = 50 x 10°, M = 0.9, and D = 2 to 22 in 4-min steps. 
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6 Test direction parameter • — * I 1 + 0.2 N?) ^ 1 * 65650 ~ C-<10 6 


Test direction A8.C, is the least energy least 'q ' locus 
Test condition set is monotonically ordered in [<p •••■£ .. 



15. Test direction design— least energy consumption locus. 
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Energy state diagram and test direction locus in cryogenic 
wind tunnels. 
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Figure 17. A cryogenic tunnel test program of 38 flow parameter conditions, 

R -M-D. 
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Figure 18. Order of execution of the cryogenic tunnel test program. 
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Figure 19. Minimum energy test direction design for the cryogenic tunnel 
test program. 
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